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ABSTRACT 

Although the widely-used analytical afterglow model of gamma-ray bursts (GRBs) predicts a sharp 
cooling break Vc in its afterglow spectrum, the GRB observations so far do not indicate a clear 
evidence for a cooling break in their spectra or its corresponding temporal break in their light curves. 
Employing a Lagrangian description of the blast wave, we conduct a sophisticated calculation of the 
afterglow emission. We precisely follow the cooling history of non-thermal electrons accelerated into 
each Lagrangian shell. We show that a detailed calculation of afterglow spectra does not in fact give 
rise to a sharp cooling break at Vc- Instead, it displays a very mild and smooth transition, which 
occurs gradually over a few orders of magnitude in energy or frequency. The main source of this slow 
transition is that different mini-shells have different evolution histories of the comoving magnetic field 
strength B, so that deriving the current value of Vc of each mini-shell requires an integration of its 
cooling rate over the time elapsed since its creation. This is in contrast to the un-justified analytical 
method, which takes the current i? as a constant value over the cooling history. We present the time 
evolution of optical and X-ray spectral indices to demonstrate the slow transition of spectral regimes, 
and discuss the implications of our result in interpreting GRB afterglow data. 

Subject headings: gamma-ray burst: general — radiation mechanisms: non-thermal — shock waves 



1. INTRODUCTION 

The broad-ba nd afterglow emis si on of gamma-ray 
bursts (GRBs) (jCosta et all 119971 : Ivan Paradijs et al.l 
[TqqI . has been interpreted as synchrotron radiation 
from a relativistic blast wav e, which sweeps up a 
surrounding am bient medium (jMeszaros fc ReesI [19971: 
iSari et ahl 119981 ). This widely- us ed analytical model of 
GRB afterglows (jSari et al.lll998[ ) predicts a sharp cool- 
ing break at Vc in their spectra, determined by syn- 
chrotron cooling rate of non-thermal electrons during the 
dynamical time scale of the relativistic blast wave. 

Although the existence of a sharp cooling break is one 
of the main ingredients of the analytical model, the GRB 
observations so far do not indicate a clear evidence for the 
cooling break in their spectra or its corresponding tem- 
poral break in light curves with a spectr al index change 
across the break. The Swift satellite ( Geh rels et al.ll20M ) 
has accumulated 8 years of afterglow data. In the X-ray 
band where Vc likely shows up, usually one or two steep- 
ening breaks are observed in the canonical X-ray light 
curves ([Zhang et al.ll2006l: iNousek et al.ll2006[) . However, 
the change of decay slope is inconsistent with Vc crossing, 
and more importantly, there is essentially no spectral in- 
dex c hange across the temporal breaks (e.g. ILiang et"al] 
[200l . The missing Vc in individual GRBs may be under- 
stood as that the shock parameters are such that there 
is no Vc crossing during the observational time, but the 
lack of Vc in the entire Swift afterglow data set suggests 
that this must be an effect intrinsic to G RB ph ysics. 

The simple analytical model (jSari et al.lll998[ ) assumes 
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that the entire postshock material forms a single zone 
with the same energy density and magnetic field. The 
single shocked zone is endowed with a broken power-law 
electron energy dist ribution with a break at a cooling 
Lorentz factor 7c. iGranot fc Saril (|2002[ ) relaxed this 
assumption and introduced other factors such as th e 
Blandford-McKee solution (jBlandford fc McKed [TQTI . 
the curvature effect, and adiabatic cooling to describe 
the shocked material, and showed that the shape of spec- 
tral breaks is not as sharp as predicted by the analytical 
model. 

In this Letter, we perform a detailed study on the for- 
mation of GRB afterglow spectra, by ado pting a La- 
grang ian description of the shocked region (lUhm et al.l 
120 121 hereafter U12). We obtain a very smooth Vc in the 
afterglow spectrum. We identify its main physical origin, 
and show that a smooth Vc is ubiquitous and intrinsic to 
a wide range of astrophysical phenomena involving syn- 
chrotron cooling. 

2. SMOOTH COOLING BREAK AND ITS PHYSICAL 
ORIGIN 
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Following lUhml (|201lD and U12 (Section 2), we make 
use of a semi-analytic formulation of a relativistic blast 
wave in order to find its dynamics accurately from an ini- 
tial coasting phase through a deceleration stage. Then, 
as described in U12 (Section 3), we adopt a Lagrangian 
description of the blast wave to conduct a sophisticated 
calculation of its afterglow emission. The blast is viewed 
as being made of many different Lagrangian shells {(5m*}. 
Here the index i is used to denote each Lagrangian shell. 
We keep track of an adiabatic evolution of each shell (5m* 
to find a time evolution of the magnetic field of the 
shell (U12, Section 3.1). The cooling history of an elec- 
tron spectrum is also followed individually for every shell 
(5m* (U12, Section 3.2). 

Radiative and adiabatic cooling of an electron with 
the Lorentz factor 7e is described by the first and second 
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term below, respectively, (U12, Equation (17)) 

A (1) ^ _ 1 (]_\ d}np^ 

dt' \7e/ GTTnieC 4 \7ey dt' 

Here t' is the time measured in the co-moving fluid frame, 
ax is the Thomson cross section, rrie is the electron mass, 
c is the speed of light, and and are respectively the 
magnetic field and the pressure of the shell 6m^ where 
the electron sits in. The Compton parameter Y that de- 
scribes a contribution of inverse Compton scattering to 
the cooling rate of electrons is omitted here. We assume 
that electrons are accelerated into a power law distribu- 
tion, and for the shell Sm^, the minimum Lorentz factor 
of the injected electrons is 7^. The differential Equation 
([1]) governs cooling of all the electrons inside the electron 
spectrum. Thus, we use Equation ([T]) to find the time 
evolution of both 7^ and the "cooling" Lorentz factor 7* 
for each Lagrangian shell (5to*, as the blast wave propa- 
gates. 

When the shell 6m^ is created at a shock front (forward 
shock (FS) or reverse shock (RS)) at the co-moving time 
t'^, the electron spectrum injected into the shell (5m* has 
the initial Lorentz factors 

llnitd^^iniit'i), 7:(i0=+OO, (2) 

where 

7inj(t')=l + ^^ee[7p(0-l]- (3) 

p~ 1 me 

Here, p is the slope of the electron spectrum, is the 
proton mass, and 7^ is the mean Lorentz factor of protons 
in the fresh shell created at the shock front (FS or RS) 
at co-moving time t'. Then, the Lorentz factors Jmi^'j) 
and 7c (ij) of the shell Sm"^ at a later time t'^ (> i-) can 
be found by integrating Equation ([T]) from t'^ to t'j . 

We take the following simple fiducial example to cal- 
culate the afterglow spectra. (1) A constant density 
pi{r)/mp = ni(r) = 1 cm~'^ is assumed for the ambient 
medium. (2) The ejecta has a constant kinetic luminos- 
ity Lcjir) — Lq — 10^^ erg/s for a duration of Tb = 5 s, 
so that the total isotropic energy of the burst is to be 
Eb — Lq Tb = 5 X 10^^ erg. (3) The ejecta is assumed to 
emerge with a constant Lorentz factor Foj = 300, so that 
the RS is short-lived. (4) The burst is assumed to be 
located at a redshift z — 1. (5) The microphysics param- 
eters are adopted as p — 2.3, €e = 10~^, and €b — 10~^ 
(for slow cooling) or 10~^ (for fast cooling). 

We find the dynamics of this blast wave by making 
use of Uhm (2011) and U12 (Section 2). Since the RS 
is short-lived here, we focus on the FS afterglow spectra 
from the FS shocked region (i.e., region 2). The magnetic 
field B'^{t') of 5to* in region 2 at any co- moving time t' 
(t'i <t' < t'j) is found following U12 (Section 3.1) while 
assuming, for simplicity, a constant pressure profile over 
region 2, i.e., p^{t') = Pf{t') for all shocked shells {5m*} 
in region 2 at every co-moving time t' . Here, Pf is the 
pressure at the FS front. Integrating Equation (jlj from 

to tj with the initial values (Equation ([2])) at t-, we 
obtain the Lorentz factors 7m(^j) and 7c(ij) of each shell 
(5to' at time t'^. 



Taking also into account the equal-arrival-time "cur- 
vature effect" of spherical shells, the Doppler boosting of 
radial bulk motion, and different radii {r' } of individual 
shells (U12, Section 3.3), we find the FS afterglow spec- 
tra for the fiducial example above. The result displays 
a very mild and smooth transition over a few orders of 
magnitude in frequency or energy in the vicinity of a 
cooling break. 

In an effort of identifying the main source of this 
smoothness, we first remove the curvature effect and 
place all the shocked shells at the same radius rf where 
the FS is located. The result is shown in Figure [H For 
panels (a) and (b), we use cb = 10~^ (slow coohng), 
and for panels (c) and (d), we use = 10^^ (fast cool- 
ing). As shown in panel (a), the FS afterglow spectra 
still exhibit a very mild and smooth cooling break. In 
particular, the spectral indices shown in panel (b) high- 
light that a transition from /? = {p—l)/2 spectral regime 
to (3 — p/2 segment is very slow and takes several orders 
of magnitude in observer time tobs- The outcome in pan- 
els (c) and (d) (fast coohng) is even more interesting. 
An expected transition from (3 — —1/3 to /? — 1/2 is 
not present. The characteristic spectral segment with 
/3 = 1/2 is not even reproduced. Instead, an almost flat 
(/3 « 0) spectral segment is observed at around the opti- 
cal bands. A more detailed investigation regarding this 
fast cooling spectra will be presented in an accompanying 
paper. 

Since the curvature effect is not the main source of 
smoothing the cooling break, the next source we suspect 
is the adiabatic cooling term included in Equation ([T]). 
Dropping out this term from Equation ^ and integrat- 
ing it from t'^ to t'p we get 

for 7*(ij) at time t'^ since ^l{t[) = +00. For the same 
fiducial example, we calculate the FS afterglow spectra 
while making use of Equation (|4]) instead of Equation 
([T]) in the calculation of 7c(ij)- The resulting afterglow 
spectra are not very different from those shown in Fig- 
ure [TJ Thus, the adiabatic cooling term in Equation ([l} 
is not the main source of this slow transition through the 
cooling break. 

The analytical method of lSari et all (|1998D uses the in- 
stantaneous B to estimate the cooling time scale, which 
implicitly assumes that the magnetic field strength B did 
not evolve during the dynamical evolution of the blast- 
wave. To check whether this is the main source of dis- 
crepancy, we take the B^ term out of the integration in 
Equation Q. Using its current value at time t^-, we get 

which resembles the widely used expression for 7c of Sari 
et al. (1998), their Equation (6). Employing Equation 
([5]) in the calculation of 7c (ij) and for simplicity, adopt- 
ing 7m (ij) = 7inj(ij) (Eq.(3)) in all shocked shells {(5m*}, 
we find the FS afterglow spectra for the fiducial example. 
The result is shown in Figure [2l Panel (a) now displays 
a considerably sharp cooling break. Thus, using Equa- 
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tion ([S|), we closely reproduce the analytical afterglow 
spectra shown in Sari et al. (1998). On the other hand, 
panels (c) and (d) indicate that the predicted spectral 
segment with /3 = 1/2 in the fast cooling regime is still 
not recovered. This will be investigated further in the 
accompanying paper. 

Thus, we have identified the main source of the smooth 
cooling break. The approach made above to get Equation 
([5|) from Equation Q cannot be justified, since the mag- 
netic field in the shell (5m* must be a time-dependent 
quantity, which requires an integration of ^ over time 
in order to correctly follow the cooling history of elec- 
trons in the shell Sm'^, as shown in Equation Q. 

3. CONCLUSIONS AND DISCUSSION 

In this Letter, we perform a detailed study on the 
formation of GRB afterglow spectra, by adopting a 
Lagrangian descr iption of the shocked region following 
lUhm et aD (|2012[ ). We precisely follow the cooling his- 
tory of the electron spectrum for each individual La- 
grangian shell, and integrate over all shells to find the 
instantaneous flux spectra. We show that this detailed 
calculation gives rise to a very mild and smooth cooling 
break, which occurs gradually over a few orders of magni- 
tude in energy. We identify the main source of this slow 
transition as due to the different B evolution histories 
of different mini shells. This gives rise to an additional 
spreading of Vc for different shells, aside from the simple 
age difference {t'j —t'^). This extra spreading is the main 
source of the smooth Vc- 

We have shown that this effect exists regardless of 
whether the curvature effect of a relativistic spherical 
shell is taken into account. It does not depend on the de- 



tails of the blastwave dynamics or whether it is FS or RS. 
In fact, it is an intrinsic effect relevant to a wide range of 
astrophysical phenomena that invokes synchrotron cool- 
ing of electrons. 

The smooth Vc presented in Figure [T] only shows the 
"minimum" smoothness. Adding the equal-arrival-time 
curvature effect or the Blandford-McKee profiles would 
introduce more spreading and can only make the break 
smoother. 

It is interesting to note that the injection break v„i 
is much sharper than Vc (although it can be smoother 
when the curvature effect is taken into account). If a 
sharp spectral break is observed in an astrophysical phe- 
nomenon, this break is very likely an injection break, and 
cannot be a cooling break. 

Simple analytical GRB afterglow models predict some 
"closure relations" between the temporal decay index a 
and the s pectral index /3 (e.g. iMeszaros fc ReesI fT"997l: 
Sari et al.iri 998. 1999; ChcvaheL&Li||2000; Dai & Chengl 
20011: iZhang fc Meszaros.,20041 ). These models have been 
applied to the Swift XRT data (from which both tempo- 
ral decay index and spectral index can be extracted) to 
test t he validity of the afterglow models (e.g.HTiang et all 
20071 '2008: Panaitescul 120071 : iWillingale et all l2007t 
Evans et al.. ,2009^) . The result of this paper suggests 
that the analytical closure relation above Vc is usually 
not achieved. The data could still be consistent with the 
afterglow theory even if the data fall into the grey zone 
between the analytical v < Vc and v > closure relation 
lines in the a — (3 plane. 

We thank KIAA for hospitality where this research is 
carried out. BZ acknowledges a Cheung Kong Scholar- 
ship in China. This work is partially supported by NSF 
AST-0908362. 
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Fig. 1. — Afterglow spectra of the fiducial example at different observation times. Equation (1) is used to calculate the cooling history of 
electron energies. The curvature effect is not taken into account to allow identification of the key physical origin of Vc smoothing. Following 
parameters are adopted: L^j{T) = Lq = 10^^ erg/s, = 5 s, = 5 X 10^^ erg, ni = 1 cm" ^, F^j = 300, z = 1, p = 2.3, £e = 0.1. Top 
Left: slow cooling case with eg = 0.001. Bottom Left: spectral index 13 evolution as a function of observer time t^^^^ for the slow cooling 
case, which shows a very slow transition of Uc crossing. Top Right: fast cooling case with eg = 0.1. Bottom Right: /3 evolution for the fast 
cooling case. 
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Fig. 2. — Calculation of the same problem as Figure [T] but with the simple, unjustified prescription as described in Equation The 
notations of all 4 panels are the same as Figure [T] A sharp Uc is reproduced. 



